Dynamical role of system-environment correlations in non-Markovian dynamics 
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We analyse the role played by system-environment correlations in the emergence of non-Markovian dynam- 
ics. By working within the framework developed in Breuer et al, Phys. Rev. Lett. 103, 210401 (2009), we 
unveil a fundamental connection between non-Markovian behaviour and dynamics of system-environment cor- 
relations. We derive an upper bound to the rate of change of the distinguishability between different states of 
the system that explicitly depends on the establishment of correlations between system and environment. We 
illustrate our results using a fully solvable spin-chain model, which allows us to gain insight on the mechanisms 
triggering non-Markovian evolution. 
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The study of non-Markovian quantum dynamics is gath- 
ering substantial interest due to key advances in the analy- 
sis, understanding and even simulation of non-trivial system- 
environment effects I id. A particularly significant step for- 
ward in this context has been performed with the formula- 
tion of new theoretical tools able to characterize and quantify 
the deviations of given dynamics from Markovianity ||2|-|6|]. 
Each of such instruments addresses a specific manifestation 
of non-Markovianity, therefore embodying iti principle a dis- 
tinct quantitative measure. Although the relationship among 
such diversified approaches has been considered |]7J| and a full 
reconciliation seems foreseeable ISO, the fundamental respon- 
sibles for the occurrence of non-Markovian features remain 
largely elusive. 

In particular the question of whether system-environment 
correlations (SECs) are of any importance in the emergence 
of non-Markovian dynamics is essentially open. Although a 
possible role of SECs was hinted in some formulation of non- 
Markovian dynamical maps llOfl, their connection with quan- 
titative measures of non-Markovianity has not been explored, 
to the best of our knowledge. Here, we aim at filling this gap 
by formulating a theory that makes such a connection explicit 
and indeed experimentally testable. We show that, by making 
use of the properties of the trace distance, it is possible to es- 
tablish a quantitative link between the non-Markovian nature 
of a process and the existence/evolution of dynamical SECs. 
Our key result is the formulation of an upper bound for the 
derivative of the trace distance between two evolving states of 
the system that depends explicitly on SECs and environmental 
distinguishability. In turn, this bound can be used to witness 
the occurrence of SECs by monitoring experimentally the be- 
havior of the trace distance. The development of increasingly 
accurate and reliable techniques for manipulation of quantum 
systems in photonics and condensed-matter physics will soon 
open up the possibility to harness the interaction between sys- 
tem and environment. In this perspective, our results provide 
an analytic tool to gather useful information on the way a sys- 
tem and its environment share quantum correlations by moni- 



toring a figure of merit of simple experimental access. 

In order to link our formal findings to an interesting phys- 
ical case, we use an exactly solvable quantum many-body 
model. We consider a quantum spin chain ruled by XX-like 
inter-particle couplings in a transverse magnetic field, a sys- 
tem that has been recently used to ascertain non-Markovian 
dynamical features llllll . This provides a pragmatic scenario 
where we illustrate our results: on one hand, we test the tight- 
ness of the bound using a physically motivated example. On 
the other hand, we quantify the amount of SECs and show the 
similarity between the behavior of the rate of change of SECs 
and that of the derivative of the trace distance. As we dis- 
cuss thoroughly, our investigation almost naturally paves the 
way to the critical assessment of the definition of "information 
flow" upon which the measure in ||4[] is built. 

Trace distance-based measure of non-Markovianity. — We 
start by setting up a framework for relating non-Markovianity 
and SECs and based on state distinguishability, which is mea- 
sured, throughout this paper, by the trace distance D{p\,p2) — 
\\P\ ~ P2II/2 between two states pi and p2- Here, ||A|| — 
TrVi?A is the trace- 1 norm of a matrix A. Given a quan- 
tum system that is prepared in state pj {j = 1,2) with 
probability 1/2, the average probability for an observer 
to guess correctly such preparation is [1 -t- D(pi,p2)]/2. 
In what follows, we exploit the subadditive property of 
the trace norm with respect to the tensor product, i.e. 
D{pi®p2,p3®p4)<D{pi,p3)-\-D(p2,P4) withpt (^=1,..,4) ar- 
bitrary density matrices, and its contractivity under a positive 
trace-preserving map ^, D{^',[pi],^',[p2]) < D(pi,p2). Con- 
tractivity gives rise to the decrement of the trace distance un- 
der completely positive (CP) maps and is a key ingredient in 
the definition of non-Markovianity. 

According to Ref. |]4[], a process is defined as non- 
Markovian if there is a pair of initial states pi,2(0) of the sys- 
tem and a time t in their evolution such that 

(r(t,pi,2m = ^D(pi{t),p2(t))>0. (1) 

dt 

The associated measure of non-Markovianity is then calcu- 



lated as 7V= maXp, ,(0) J" dt cr(t,p 1^2(0)), where a+ is the 
union of the temporal domains where cr(f,pi 2(0))>0. The 
optimization over the pair of input states ensures that such a 
figure of merit faithfully reveals non-Markovianity. 

Bound to non-Markovianity.- We now provide the con- 
nection between the emergence of non-Markovianity as wit- 
nessed by Eq. (HJ and the evolution of SECs by formulating 
an upper bound to the derivative of the trace distance between 
two system states p^ and p^ at a generic instant of time t. Our 
result can be effectively formulated as the following 
Theorem. For any quantum process described by a com- 
pletely positive map with associated system-environment in- 
teraction ruled by the propagator U(t,0) — e^'^' we have 



-(^Pi,2)^5(min||Tr£{[i/,pf(0®(pf(r)-pf(r))]}|| 
+ \\rrE{[H,(xl^(t)-//(t)) 



where 



Pj2 cifs arbitrary systems states, X iO 
are the SECs, 



(2) 



and 



.J ,., p^(f) ® p;(0 

pf^>(0=Tr£(5)[f/(,,,opf(f')f/J,,,)] (for any t'<t) is the 
system (environment) state at the time t. 

The technical steps needed to prove this statement are 
sketched in the Appendix. Here we focus on the signifi- 
cance and the implications of this Theorem, which is the cen- 
tral result of our study. Eq. (|2) bridges the occurrence of 
non-Markovianity (as revealed by a growing trace distance) 
with the dynamics of SECs and environmental distinguisha- 
bility. The upper bound consists of two different contribu- 
tions. The first term contains the information about the two 
different states of the environment. By decomposing the in- 
teraction Hamiltonian in terms of the eigen-operators A^ and 
Bf , so that H= 2ff ^a®^ff ' '^he first term can be recast into the 
form II [Za7aAa,pl(i)(t)] II with y„= + Tr£{B„(pf (f)-pf (f))). 
Here y clearly depends on the difference between the environ- 
ment states. The second term of the bound accounts for the 
presence and evolution of SECs. This term contains all the 
non-diagonal elements of the S -E state, depends by definition 
on both the reduced state of the system and the total system- 
environment state, and is such that Tis(E)[x^-^{t')] = 0. 

The bound in Eq. ^ shows that if SECs are not produced 
across the evolution and the environment in left in the same 
conditions regardless of the state of the system, the process is 
necessarily Markovian. This demonstrates the intimate con- 
nection between SECs and the changes in distinguishability 
of different input states in a non-Markovian evolution, broad- 
ening the view on the occurrence of such effects. On the other 
hand, our upper bound states that the creation of correlations 
can be compatible with a Markovian dynamics. 

An important point to stress is that the Theorem is formu- 
lated under the assumptions of zero initial SECs; the ini- 
tial state of system and environment is factorized. A unifi- 
cation between witness of initial SECs and witness of non- 
Markovian dynamics is proposed and thoroughly discussed in 
Ref. H. 



Physical Model. — ^We now test Eq. ^ against a physically 
motivated situation that will help us illustrate its deep phys- 
ical implications. We consider the generation of SECs and 
their relation with non-Markovianity in a simple unidimen- 
sional quantum many-body system embodied by A^ -H 1 spin- 
1/2 particles (labelled n - 0, l,..,N) mutually coupled via 
an XX model and subjected to a transverse magnetic field. 
We are interested in the bipartition consisting of the small 
system given by spin and the environment represented by 
the remaining A^ particles. Assuming units such that h - I 
across the manuscript, the corresponding Hamiltonian model 
is H-HsE + He with 



HsE = -2yo(o"oO"i + o-qO-i) 



N 



He = -27 y {&'„&:,, + d-ld-l^,) -IBVcji 



(3) 



n=l 



«=1 



where o-f, is the A:-Pauli matrix {k=x,y,z) for particle n, B 
is the amplitude of the magnetic field affecting S and J 
(7o) is the inter-environment (system-environment) coupling 
strength. This model presupposes that the free evolutions of 5 
and E are identical, thus allowing the passage to the interac- 
tion picture without the introduction of time-dependent coeffi- 
cients. The open-system evolution of S and its consequences 
for non-Markovianity have been analysed in Ref. llllll . where 
it was found that, for interaction times that are well within the 
recurrence time of the system (when any information prop- 
agating across the chain returns to S after reaching the end 
of the chain), there is a working point defined by (Jq/J, B/J) 
at which the measure of non-Markovianity N is null. As 
the optimization inherent in the definition of such measure 
is achieved for system states lying on the equatorial plane of 



the Bloch sphere llllll . we consider the input states p^2^'^) ~ 



|±)<+| with |+)=(|0) + |1))/ V2. The environment is initial- 
ized in pf 2(0) = p^„. - ®J^7' |0),<0|, so that the total system- 
environment state is p^f(O) - pf 9(0) ®p^;. 

Fig. [r|(a) compares cr(t, pi,2) and its upper bound under the 
above initial conditions and for Jq/J-I, B/J - 0.01. We 
have restricted the width of the time-window that we consider 
to values within the time at which finite-size effects are ex- 
pected to occur. This is easily estimated by considering that 
the maximum single-excitation group velocity of our model 
is 2, which implies that it takes at least a time A^ for a sin- 
gle excitation to leave S and come back to it. As the right 
hand side of Eq. ^ is non negative, a quantitative compari- 
son between the derivative of the trace distance and its upper 
boundismeaningfulonly when cr(f,pi 2(0) > 0. Interestingly, 
the bound becomes very tight as soon as the derivative of the 
trace distance becomes positive. 

The two terms appearing in the bound have a different time 
behavior: the one depending on the reduced environmental 
states contributes the most to the right hand side. The second 
term, which is associated with the correlations, has local min- 
ima whenever cr(f,pi 2) crosses the horizontal axis. This ob- 
servation prompts us to study the behavior of other quantities. 
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FIG. 1: (Color online) (a) We show cr(f,pi_2) (solid blue line) and 
the right hand side of Eq. ([2j (which we label S for convenience 
and display as a dashed orange line) against the rescaled interaction 
time Jt for a chain of A' = 9 particles (i.e. an environment with 8 
spins) when 5 is prepared in either |+) or |-) (see inset for a sketch 
of the physical situation at hand), (b) We compare the behaviour 
of cr(t,pi2) (solid blue line) with the evolution of the environmental 
indistinguishabUity £= 1 -D(p^(i),p'^{t)) (amplified by a factor 10 for 
easiness of visualization, dotted red line) and the quantity X defined 
in the body of the manuscript (dashed green line). In both panels, 
Jq/J=1 and B = IO""y, while shaded regions highlight the regions 
of non-Markovianity. 



In Fig.[T|(b) we plot the evolution of the trace distance of the 
correlations X=D(x^^^ (t), x^^ (t)) and the indistinguishabUity 
between the environmental states S-l-D(p^{t),p^(t)) find- 
ing that they are both minimal when cr(f,pi 2) changes sign 
from being negative, thus showing that non-Markovianity oc- 
curs when the environmental states become perfectly distin- 
guishable [i.e. D(p^(f),pf(0)=l]- At these times, a careful 
analysis shows that not only X - but also the correlations 
within each state almost disappear (|[;\f'^f (Oil - 0), thus leaving 
system-environment states that are basically factorized. 

Time-derivative of quantum mutual information. — We 
know pass to the quantification of SECs in non-Markovian 
dynamics. The quantum mutual information (I) between 
S and E quantifies the total amount of information shared 
between system and environment, and accounts for classi- 
cal and quantum correlations between the parts. Given that 
the state of the total system (plus environment) is pure and 
evolves via unitary dynamics, the quantum mutual infor- 
mation is simply equal to twice the von Neumann entropy 

>S,;v(0=-Tr[p'^,(01og2P'^2(^)] of the system (coinciding also 
with the entanglement between S and E). Notice that the 
von Neumann entropy is the same for the two system's states. 
Fig. [2] (a) compares the derivative of the von Neumann en- 
tropy with cr(r,pi 2): with the exception of the initial part 
of the dynamics, the two functions exhibit the same qualita- 
tive behaviour When the evolution starts the quantum mu- 
tual information increases, so SECs are created, leading to 
a loss of information over the system's state and thus to de- 
coherence. Later on, both mutual information and SECs de- 
crease: this is the time-window where re-coherence takes 
place. If such a process becomes significant, and thus cor- 
relations become small enough (implying an actual comeback 
of the information that was previously encoded in the S -E 
state to S only), the system "jumps" to the non-Markovian 
regime with an abrupt change of cr(f,pi 2). Quite interest- 
ingly, such re-coherence processes do not coincide with the 




FIG. 2: (Color online) (a) We show cr{t,p[2) (solid blue line) and 
■SyN(t) (dashed red curve) against the interaction time Jt for Jq/J = 
1,6/y = 10"-. The trace distance D(p^^(t),p2(t)) and entanglement 
associated to this situation are reported for the same time- window. 
The shaded regions highlight non-Markovianity. (b) Same analysis 
but for the Markovian point corresponding to Jo/J=l and B/J=l/2. 



non-Markovianity period, but seem to act as precursors of it. 

This behavior is typical of the working points at which the 
model at hand gives rise to non-Markovian processes certi- 
fied by Ni^Q. On the other hand, as mentioned above, a 
parameter-regime exists for which the many-body dynamics 
corresponds to a perfectly Markovian process for spin S Jll l. 
This is reported in Fig. |2](b). Clearly, the much higher de- 
gree of quantum correlations between S and E at the times 
at which, in panel (a), the non-Markovian thresholds were 
passed, leaves cr(f,pi 2)<0. The quantum mutual information 
approaches zero (and cr(f,pi 2) is positive again) only for time- 
windows wide enough to cover the effects of to the chain's fi- 
nite size. However, this sort of non-Markovianity should be 
set apart from the features discussed here and associated to 
shorter timescales. While the former are due to the physical 
return of excitations to the S particle, the latter have a much 
deeper origin, now captured by Eq. ([2) and our analysis. 

Conclusive Remarks. — We have proposed and analysed in 
depth a new approach to the establishment of non-Markovian 
dynamics, based on the assessment of the state of the envi- 
ronment and the explicit presence of SECs during the open 
evolution. Our perspective is that, in order to get insight into 
the reasons behind a non-Markovian process, one should fo- 
cus on the structural modifications that the interaction with S 
induces on the system-environment states. Quite significantly, 
our investigation appears to support such a vision by linking 
explicitly one of the measures for non-Markovianity to the 
presence of SECs and possible deviations of the environment 
from its initial state. We have provided an upper bound to 
the rate of change of the distinguishability between different 
states of S that clearly shows a dependence on such key fea- 
tures. The tightness of the bound was investigated using a 
spin-chain system whose rich non-Markovianity diagram al- 
lowed the exploration of various interesting dynamical condi- 
tions lull and the study of the subtle role played by SECs in 
the emergence of non-Markovianity. 

We believe that SECs can be instrumental also to ad- 
dress critically the interpretation of the measure of non- 
Markovianity based on the trace distance. Such a quantifier is 
operationally interpreted in terms of the so-called "flow of in- 



formation" to and from the environment: for non-Markovian 
systems, the natural direction of the decoherence process lead- 
ing the system to lose energy/coherence/quantumness into the 
environment is, in certain time windows, reversed. The re- 
versal of the flow of information is assumed to happen when 
the derivative of the trace distance increases. Nevertheless, 
when used in the context of non-Markovian dynamics, the 
term information-flow does not have a uniquely defined math- 
ematical meaning. Our view is that SECs embodies a proper 
figure of merit to define information flow. In such a perspec- 
tive, information would be identified by the quantum mutual 
information I shared by system and environment, thus quan- 
tifying the total amount of SECs. On the other hand, its flow 
will be encompassed by the time derivative of I . An inflow 
(outflow) of information to (from) the system would be sig- 
naled by its negative (positive) time derivative. The applica- 
tion of such a definition to the physical model at hand would 
challange the commonly accredited view according to which 
S recovers information while the trace distance is positive. In 
fact, in the physical model studied here, the re-coherence pro- 
cess, which is associated with an increased purity of the sys- 
tem's state and the decrease of SECs, occurs before the tran- 
sition to a non-Markovian window and, in a sense, precedes 
it. 

Our study calls out loud for a clearer definition of infor- 
mation flow in open quantum systems and proposes the us 
of correlations between system and environment as the figure 
of merit against which such an analysis should be performed. 
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Proof of the stated Theorem 

We now provide the steps needed in order to derive the up- 
per bound to cr(f,pi 2) stated in Eq. (|2]i. As discussed in the 
main body of the paper, we consider a quantum process de- 
scribed by a completely positive map and take the derivative 
of the trace distance between two system states, p^ and p^ , at 
a generic instant of time t. We get 

1 ||pf(f)-p^(f)||-||pf(0-pf(OII ^,, 

o-(f,pi,2) = - lim — ! '- -— -! ^ . (4) 

2 I'^i t - f 

We now use the decomposition of the S -E state at time t' 
p^.^(f) = p^. it') ®p^{t') +x^-^{t') and the triangular inequality 
for the trace distance to obtain 

„S/^x „S/.Mi , i|rr„ rA /„Sfj 



\\p\(t)-pim < \\TxE[U(,fM^t')-pW)) ®pt(2,(f')c/(;,,,)) 



+ min \\TxE{U(,,n{pl{t')) ® (pf (f') - pWwI,,) 

+ IITr£{?/(,,oOrf (f') -xt^(t'))ul„^}\\- 



(5) 



By using the fact that the evolution of a system through the 
unitary operator U(,fi) is described by the completely posi- 
tive (CP) dynamical map, the first term in the right hand side 
of Eq. © is written as ||0(,,,')[p^(/')] - 3)(,_,')[p|(f')]||, where 

^(,,nlPi(t')'i = Tr£{f/(,,,)pf (f) ® pf(2)(f')f/(/,/')) is a CP trace- 
preserving map. Using again contractivity we have that 



mrrMif')] - <^(r,r'M(f')]\\ < \\p\(t') - p^t' 



(6) 



By plugging Eq. ^ into Eq. (|5]l and substituting back in 
Eq. (|4]i, we find 



0-{t,pi,2) ^ 



1 , (I|1>(,,r)[pf(f')-P^(0]l|-Ilpf(0-P^(0 



ilimli 



t - f 



l|Tr£{f/(,,oC^f(f')-;rf(f'))<,,) 



t-f 



+ min 

k=l,2 






(7) 

The first term in the right hand side is non-positive due to 
the contraction property of the trace distance under CP maps . 
We discard it, thus providing a loser upper bound to cr(f,pi 2). 
In order to manage the second and third terms we assume to 
know the Hamiltonian H regulating the S -E interaction and 
expand U(t.t') - e"'^*'"'' in power series, stopping it at first 



order in H. We find 



lim 



\\T,E{Ui,,n(j(\Ht')-xlHt'))Ul,^ 

(t-n 



An analogous calculation holds for the third term, from 
which Eq. ^ is then straightforwardly found. 



